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ABSTRACT 

Recent time-dependent, ideal-magnetohydrodynamic (ideal-MHD) simula- 
tions of polar magnetic burial in accreting neutron stars have demonstrated that 
stable, magnetically confined mountains form at the magnetic poles, emitting 
gravitational waves at /* (stellar spin frequency) and 2/*. Global MHD oscilla- 
tions of the mountain, whether natural or stochastically driven, act to modulate 
the gravitational wave signal, creating broad sidebands (full-width half-maximum 
~ 0.2/*) in the frequency spectrum around /* and 2/*. The oscillations can en- 
hance the signal-to-noise ratio achieved by a long-baseline interferometer with 
coherent matched filtering by up to 15 per cent, depending on where /* lies 
relative to the noise curve minimum. Coherent, multi-detector searches for con- 
tinuous waves from nonaxisymmetric pulsars should be tailored accordingly. 

Subject headings: gravitation — gravitational waves — stars: magnetic fields 
— stars: neutron — stars: rotation 

1. Introduction 

Nonaxisymmetric mountains on accreting neutron stars with millisecond spin periods 
are promising gravitational wave sources for long-baseline interferometers like the Laser 
Interferometer Gravitational Wave Observatory (LIGO). Such sources can be detected by 
coherent matched filtering without a computationally expensive hierarchical Fourier search 
(Brady et al. 1998), as they emit continuously at periods and sky positions that are known 
a priori from X-ray timing, at least in principle. Nonaxisymmetric mountains have been 
invoked to explain why the spin frequencies /* of accreting millisecond pulsars, measured from 
X-ray pulses and/or thermonuclear burst oscillations (Chakrabarty et al. 2003; Wijnands 
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et al. 2003), have a distribution that cuts off sharply above /* 0.7 kHz. This is well below 
the centrifugal break-up frequency for most nuclear equations of state (Cook et al. 1994), 
suggesting that a gravitational wave torque balances the accretion torque, provided that the 
stellar ellipticity satisfies e ~ 10 -8 (Bildsten 1998). Already, the S2 science run on LIGO 
I has set upper limits on e for 28 isolated radio pulsars, reaching as low as e < 4.5 x 10 -6 
for J2124— 3358, following a coherent, multi-detector search synchronized to radio timing 
ephemerides (The LIGO Scientific Collaboration: B. Abbott et al. 2004a). Temperature 
gradients (Bildsten 1998; Ushomirsky et al. 2000), large toroidal magnetic fields in the stellar 
interior (Cutler 2002), and polar magnetic burial, in which accreted material accumulates in 
a polar mountain confined by the compressed, equatorial magnetic field (Melatos & Phinney 
2001; Payne & Melatos 2004; Melatos & Payne 2005), have been invoked to account for 
ellipticities as large as e ~ 10~ 8 . The latter mechanism is the focus of this paper. 

A magnetically confined mountain is not disrupted by ideal-magnetohydrodynamic 
(ideal-MHD) instabilities, like the Parker instability, despite the stressed configuration of 
the magnetic field (Payne & Melatos 2005). However, magnetospheric disturbances (driven 
by accretion rate fluctuations) and magnetic footpoint motions (driven by stellar tremors) 
induce the mountain to oscillate around its equilibrium position (Melatos & Payne 2005). 
In this paper, we calculate the Fourier spectrum of the gravitational radiation emitted by 
the oscillating mountain. In §2, we compute e as a function of time by simulating the global 
oscillation of the mountain numerically with the ideal-MHD code ZEUS-3D. In §3, we calcu- 
late the gravitational wave spectrum as a function of wave polarization and accreted mass. 
The signal-to-noise ratio (SNR) in the LIGO I and II interferometers is predicted in §4 as 
a function of M a , for situations where the mountain does and does not oscillate, and for 
individual and multiple sources. 



2. Magnetically confined mountain 

2.1. Grad-Shafranov equilibria 

During magnetic burial, material accreting onto a neutron star accumulates in a column 
at the magnetic polar cap, until the hydrostatic pressure at the base of the column over- 
comes the magnetic tension and the column spreads equatorward, compressing the frozen-in 
magnetic field into an equatorial magnetic belt or 'tutu' (Melatos & Phinney 2001; Payne & 
Melatos 2004). Figure 1 illustrates the equilibrium achieved for M a = 1O _5 M , where M a is 
the total accreted mass. As M a increases, the equatorial magnetic belt is compressed further 
while maintaining its overall shape. 
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In the steady state, the equations of ideal MHD reduce to the force balance equation 
(CGS units) 



where B, p, p — c^p, and $(r) = GM*r/Rl denote the magnetic field, fluid density, pressure, 
and gravitational potential respectively, c s is the isothermal sound speed, M* is the mass 
of the star, and is the stellar radius. In spherical polar coordinates (r,8,<f)), for an 
axisymmetric field B = V^(r, 6)/{r sin 6) x e^, (1) reduces to the Grad-Shafranov equation 



where A 2 is the spherical polar Grad-Shafranov operator, F(ip) is an arbitrary function of 
the magnetic flux ip, and we set $o — ${R*)- In this paper, as in Payne & Melatos (2004), 
we fix F{ip) uniquely by connecting the initial and final states via the integral form of the 
flux-freezing condition, viz. 



where C is any magnetic field line, and the mass-flux distribution is chosen to be of the 
form dM/dip oc exp(— where ip a is the polar flux, to mimic magnetospheric accretion 
(matter funneled onto the pole). We also assume north-south symmetry and adopt the 
boundary conditions ip = dipole at r = i?* (line tying), ip = at 6 = 0, and dip/dr = at 
large r. Equations (2) and (3) are solved numerically using an iterative relaxation scheme 
and analytically by Green functions, yielding equilibria like the one in Figure 1. 



Vp + pV$ - (47r) _i (V x B) x B = 0, 



(1) 



AV = F'(^)exp[-($-$o)/c 2 ], 



(2) 




(3) 
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Fig. 1. — Equilibrium magnetic field lines (solid curves) and density contours (dashed curves) for 
M a = 10~ 5 M Q and tp & = O.1-0*- Altitude is marked on the axes (log scale). [From Payne & Melatos 
(2004).] 
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Fig. 2.— Normalized ellipticity e(f)/e for M„/M c = 0.16,0.80,1.6, with e = 8.0 x 1.2 x 
10~ fi , 1.3 x 10~ fi respectively for b — 10. Time is measured in uidts of the Alfven crossing time, ta. 
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Fig. 3. — (Top) Fourier transforms of the wave strain polarization amplitudes h + (f) (left) and 
h x (f) (right) for M a /M c = 0.16 (dashed) and 0.8 (solid), compared with the LIGO I and II 
noise curves h 3 / yT (see §4) (dotted). The signals for M a /M c = 0.16 and 0.8 yield SNR = 2.9 and 
4.4 respectively after 10 7 s. (Bottom). Zoomed-in view after reducing h +jX (f*) an d h +jX (2f*) 
artificially by 90 per cent to bring out the sidebands. 'S' and 'A' label the signals induced by 
sound- and Alfven-wave wobbles respectively. All curves are for a = tt/3, i = 7r/3, ip*/ipa, = 10, 
and d = 10 kpc. 
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2.2. Global MHD oscillations 

The magnetic mountain is hydromagnetically stable, even though the confining magnetic 
field is heavily distorted. Numerical simulations using ZEUS-3D, a multipurpose, time- 
dependent, ideal-MHD code for astrophysical fluid dynamics which uses staggered-mesh 
finite differencing and operator splitting in three dimensions (Stone & Norman 1992), show 
that the equilibria from §2.1 are not disrupted by growing Parker or interchange modes over 
a wide range of accreted mass (1CT 7 M < M a < 1O~ 3 M ) and intervals as long as 10 4 Alfven 
crossing times (Payne & Melatos 2005). 

The numerical experiments leading to this conclusion are performed by loading the 
output (p and B) of the Grad-Shafranov code described in Payne & Melatos (2004) into 
ZEUS-3D, with the time-step determined by the Courant condition satisfied by the fast 
magnetosonic mode. The code was verified (Payne & Melatos 2005) by reproducing the 
classical Parker instability of a plane-parallel magnetic field (Mouschovias 1974) and the 
analytic profile of a static, spherical, isothermal atmosphere. Coordinates are rescaled in 
ZEUS-3D to handle the disparate radial {c^Rl/ 'GM m ) and latitudinal (R*) length scales. The 
stability is confirmed by plotting the kinetic, gravitational potential, and magnetic energies 
as functions of time and observing that the total energy decreases back to its equilibrium 
value monotonically i.e. the Grad-Shafranov equilibria are (local) energy minima. Note that 
increasing p uniformly (e.g. five-fold) does lead to a transient Parker instability (localized 
near the pole) in which < 1% of the magnetic flux in the 'tutu' escapes through the outer 
boundary, leaving the magnetic dipole and mass ellipticity essentially unaffected. 

Although the mountain is stable, it does wobble when perturbed, as sound and Alfven 
waves propagate through it (Payne & Melatos 2005). Consequently, the ellipticity e of the 
star oscillates about its mean value e. The frequency of the oscillation decreases with M a , as 
described below. The mean value e increases with M a up to a critical mass M c and increases 
with ipa,/ip*, as described in §3.1. 

In ideal MHD, there is no dissipation and the oscillations persist for a long time even if 
undriven, decaying on the Alfven radiation time-scale (which is much longer than our longest 
simulation run). In reality, the oscillations are also damped by ohmic dissipation, which is 
mimicked (imprecisely) by grid-related losses in our work. 

To investigate the oscillations quantitatively, we load slightly perturbed versions of the 
Grad-Shafranov equilibria in §2.1 into ZEUS-3D and calculate e as a function of time t. 
Figure 2 shows the results of these numerical experiments. Grad-Shafranov equilibria are 
difficult to compute directly from (2) and (3) for M a > 1.6M C , because the magnetic topology 
changes and bubbles form, so instead we employ a bootstrapping algorithm in ZEUS-3D 
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(Payne & Melatos 2005), whereby mass is added quasistatically through the outer boundary 
and the magnetic field at the outer boundary is freed to allow the mountain to equilibrate. 
The experiment is performed for r /-R* = c^R^/GM* = 2 x 10~ 2 (to make it computationally 
tractable) and is then scaled up to neutron star parameters (r /-R* = 5 x 10~ 5 ) according 
to e oc (i?*/r ) 2 and r A oc R*/r , where r A is the Alfven crossing time over the hydrostatic 
scale height r (Payne & Melatos 2005). 

The long-period wiggles in Figure 2 represent an Alfven mode with phase speed i> A oc 
M~ 1/2 ; their period roughly triples from 100r A for M a /M c = 0.16 to 300r A for M a /M c = 1.6. 
Superposed is a shorter-period sound mode, whose phase speed c s is fixed for all M a . Its 
amplitude is smaller than the Alfven mode; it appears in all three curves in Figure 2 as a 
series of small kinks for t < 50r A , and is plainly seen at all t for M a /M c = 0.8. As M a 
increases, the amplitude of the Alfven component at frequency f A ~ 17(M a /M c ) -1 / 2 Hz is 
enhanced. By contrast, the sound mode stays fixed at a frequency /s ~ 0.4 kHz, while its 
amplitude peaks at M a ~ M c (Payne & Melatos 2005). 



In this section, we predict the frequency spectrum of the gravitational-wave signal emit- 
ted by freely oscillating and stochastically perturbed magnetic mountains in the standard 
orthogonal polarizations. 



The metric perturbation for a biaxial rotator can be written in the transverse-traceless 
gauge as hj? — h+e? H + h x e*-, where e% and e*- are the basis tensors for the + and x 
polarizations and the wave strains h + and h x are given by (Zimmermann & Szedenits 1979; 
Bonazzola & Gourgoulhon 1996) 



3. Frequency spectrum of the gravitational radiation 



3.1. Polarization amplitudes 



h 



h 



■+ 



X 



ho sin a [cos a sin i cos i cos(fit) 

— sin a(l + cos 2 i) cos(2f2t)] , 

ho sin a [cos a sin i sin (fit) 

— 2 sin a cosi sin(2fit)] , 



(4) 



(5) 
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with 1 

ho = 2GI zz en 2 /dc 4 . (6) 

Here, Q = 27r/* is the stellar angular velocity, i is the angle between the rotation axis e z and 
the line of sight, a is the angle between e z and the magnetic axis of symmetry, and d is the 
distance of the source from Earth. 

The ellipticity is given by e = \I ZZ — I yy \/I , where iy denotes the moment-of-inertia 
tensor and I = ^M*B%. In general, e is a function of t; it oscillates about a mean value e, as 
in Figure 2. Interestingly, the oscillation frequency can approach Q for certain values of M a 
(see §3.2), affecting the detectability of the source and complicating the design of matched 
filters. The mean ellipticity is given by (Melatos & Payne 2005) 

U VlO- 2 Mj \10™Gj a C 2M, V 2b J V 8M C J 

(7) 

where M c = GM :¥ B 2 R 2 J(84) is the critical mass beyond which the accreted matter bends the 
field lines appreciably, b = is the hemispheric to polar flux ratio, and is the polar 

magnetic field strength prior to accretion. For = 10 6 cm, c s = 10 8 cm s -1 , and = 10 12 
G, we find M c = 1.2 x 1O~ 4 M . The maximum ellipticity, e -> 20M C /(M*6 2 ) ~ 10- 5 (6/10)" 2 
as M a — > oo, greatly exceeds previous estimates, e.g. e ~ 10 -12 for an undistorted dipole 
(Katz 1989; Bonazzola & Gourgoulhon 1996), due to the heightened Maxwell stress exerted 
by the compressed equatorial magnetic belt. Note that e(t) is computed using ZEUS-3D for 
b = 3 (to minimize numerical errors) and scaled to larger b using (7). 



3.2. Natural oscillations 

We begin by studying the undamped, undriven oscillations of the magnetic mountain 
when it is "plucked", e.g. when a perturbation is introduced via numerical errors when 
the equilibrium is translated from the Grad-Shafranov grid to the ZEUS-3D grid (Payne & 
Melatos 2005). We calculate h x (t) and h + (t) for /* = 0.6 kHz from (4) and (5) and display 
the Fourier transforms h x (f) and h + (f) in Figure 3 for two values of M a . The lower two 
panels provide an enlarged view of the spectrum around the peaks; the amplitudes at /* and 
2/* are divided by ten to help bring out the sidebands. 

In the enlarged panels, we see that the principal carrier frequencies f — f*, 2/* are 
flanked by two lower frequency peaks arising from the Alfven mode of the oscillating moun- 



1 Our ho is half that given by Eq. (22) of Bonazzola & Gourgoulhon (1996). 
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tain (the rightmost of which is labelled 'A'). Also, there is a peak (labelled 'S') displaced 
by Af ~ 0.4 kHz from the principal carriers which arises from the sound mode; it is 
clearly visible for M a /M c = 0.8 and present, albeit imperceptible without magnification, 
for M a /M c = 0.16. Moreover, e diminishes gradually over many ta (e.g. in Figure 2, for 
M a /M c = 0.16, e drifts from 1.02e to 0.99e over 500ta), causing the peaks at / = /*, 2/* to 
broaden. As M a increases, this broadening increases; the frequency of the Alfven component 
scales as /a oc M a _1 ' 2 and its amplitude increases oc M^ 2 (see §2.2); and the frequency 
of the sound mode stays fixed at / s ~ 0.4 kHz (Payne & Melatos 2005). Note that these 
frequencies must be scaled to convert from the numerical model (r /-R* = 2 x 10~ 2 ) to a 
realistic star (r /R* = 5 x 10~ 5 ); it takes proportionally longer for the signal to cross the 
mountain (Payne & Melatos 2005). 



3.3. Stochastically driven oscillations 

We now study the response of the mountain to a more complex initial perturbation. 
In reality, oscillations may be excited stochastically by incoming blobs of accreted matter 
(Wynn & King 1995) or starquakes that perturb the magnetic footpoints (Link et al. 1998). 
To test this, we perturb the Grad-Shafranov equilibrium ^ G s with a truncated series of 
spatial modes such that 

i> = ^gs{1 + S n 5 n sin[m7r(r - -R*)/(r max - R*)] sin(m6>)} (8) 

at t — 0, with mode amplitudes scaling according to a power law 5 n = 0.25m _1 , m — 2n+l, 
< n < 3, as one might expect for a noisy process. We place the outer grid boundary at 
r mm = R* + 10ro- Figure 4 compares the resulting spectrum to that of the free oscillations 
in §3.2 for M a /M c = 0.8. The stochastic oscillations increase the overall signal strength at 
and away from the carrier frequencies /* and 2/*. The emitted power also spreads further 
in frequency, with the full-width half-maximum of the principal carrier peaks measuring 
Af m 0.25 kHz (c.f. Af 0.2 kHz in Figure 3). However, the overall shape of the spectrum 
remains unchanged. The Alfven and sound peaks are partially washed out by the stochastic 
noise but remain perceptible upon magnification. The signal remains above the LIGO II 
noise curves in Figure 4; in fact, its detectability can (surprisingly) be enhanced, as we show 
below in §4. 
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Fig. 4. — (Top) Fourier transforms of the wave strain polarization amplitudes h + (f) (left) and 
h x (f) (right) for M a /M c = 0.8 with stochastic (dashed) and natural (solid), oscillations, compared 
with the LIGO I and II noise curves h 3 / yT (see §4) (dotted) corresponding to 99% confidence after 
10 7 s. (Bottom) Zoomed in view with /i+, x (/*) and h +jX (2/*) artificially reduced by 90 per cent to 
bring out the sidebands. 'S' and 'A' label the signal induced by sound- and Alfven-wave wobbles 
respectively. All curves are for a = tt/3, i = 7r/3, ?p*/'>Pa = 10, and d = 10 kpc. 
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4. Signal-to-noise ratio 

In this section, we investigate how oscillations of the mountain affect the SNR of such 
sources, and how the SNR varies with M a . In doing so, we generalize expressions for the 
SNR and characteristic wave strain h c in the literature to apply to nonaxisymmetric neutron 
stars oriented with arbitrary a and i. 



4.1. Individual versus multiple sources 

The signal received at Earth from an individual source can be written as h(t) = 
F+(t)h+(t)+F x (t)h x (t), where F + and F x are detector beam-pattern functions (0 < |i*+,x| < 
1) which depend on the sky position of the source as well as a and i (Thorne 1987). The 
squared SNR is then (Creighton 2003) 2 



N 2 



f 

Jo 



df 



\Kf)\ 2 

s h (f) 



(9) 



where S h (f) = \h 3 / yr (f)\ 2 is the one-sided spectral density sensitivity function of the detector 
(Figures 3 and 4), corresponding to the weakest source detectable with 99 per cent confidence 
in 10 7 s of integration time, if the frequency and phase of the signal at the detector are known 
in advance (Brady et al. 1998). 

A characteristic amplitude h c and frequency f c can also be defined in the context of 
periodic sources. For an individual source, where we know a, i, F + and F x in principle, the 
definitions take the form 



fc 



f 

Jo 



df 



\h(f)\ 2 

s h (f) 



-1 -1 









dff 



\Hfjf 

Shif) 



and 



^c = |[^(/ c )] 1 / 2 . 



(10) 



(11) 



These definitions are valid not only in the special case of an individual source with a = it/ 2 
(emission at 2/* only) but also more generally for arbitrary a (emission at /* and 2/*). Using 
(4), (5), (10) and (11), and assuming for the moment that e is constant (i.e. the mountain 
does not oscillate), we obtain 



f c = UxAi + 2A 2 )/( x A 1 + A 2 ), 



(12) 



2 This is twice the SNR defined in Eq. (29) of Thorne (1987). 
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- = h [S h (2fJ]- 1 ' 2 ( X A 1 + A 2 ) 1 ' 2 mna (13) 

with A x = cos 2 a sin 2 i(F+ cos i+F x ) 2 , A 2 = sin 2 a[F+(l+cos 2 i)+2F x cost} 2 , x = S h (2f*)/S h (f*), 
and i] = Sh{f c )/Sh{f*)- In the frequency range 0.2 < / < 3 kHz, the LIGO II noise 
curve is fitted well by h 3/yt (f) = 10- 26 (//0.6kHz) Hz" 1 / 2 (Brady et al. 1998), implying 
X = 4. As an example, for (a,i) = (ir/3,ir/3), we obtain f c = 1.67/*, h c = 1.22/i and 
S/N = 2.78(/*/0.6kHz)(e/10- 6 )(rf/10kpc)- 1 . In the absence of specific knowledge of the 
source position, we take F x = F + = \j\fh (for motivation, see below). 

If the sky position and orientation of individual sources are unknown, it is sometimes 
useful to calculate the orientation- and polarization-averaged amplitude h c and frequency 
f c . To do this, one cannot assume a = ir/2, as many authors do (Thorne 1987; Bild- 
sten 1998; Brady et al. 1998); sources in an ensemble generally emit at /* and 2/*. In- 
stead, we replace \h(f)\ 2 by (\h(f)\ 2 ) in (9), (10) and (11), defining the average as (Q) = 
j j Q Qd(cosa) d(cosi). This definition is not biased towards sources with small a; we pre- 
fer it to the average {Q)2 — tt -1 Jq Jq Q dad(cosi), introduced in Eq. (87) of Jaranowski 
et al. (1998). Therefore, given an ensemble of neutron stars with mountains which are not 
oscillating, we take (F 2 ) = (F 2 ) = 1/5 and (F + F x ) = [Eq. (110) of Thorne (1987), c.f. 
Bonazzola & Gourgoulhon (1996); Jaranowski et al. (1998)], average over a and i to get 
(A x sin 2 a) = 8/75 and (A 2 sin 2 a) = 128/75, and hence arrive at f c = 1.80/*, h c = 1.31h 
and (S 2 /N 2 y/ 2 = 2.78(/*/0.6kHz)(e/10- 6 )(rf/10kpc)- 1 . This ensemble-averaged SNR is 
similar to the non-averaged value for (a,i) = (7r/3, 7r/3), a coincidence of the particular 
choice. 

Our predicted SNR, averaged rigorously over a and i as above, is (2/3) 1//2 times smaller 
than it would be for a = n/2, because the (real) extra power at /* does not make up for the 
(artificial) extra power that comes from assuming that all sources are maximal (a = n/2) 
emitters. Our value of h c is 9/10 of the value of h c quoted widely in the literature (Thorne 
1987; Bildsten 1998; Brady et al. 1998). The latter authors, among others, assume a = tt/2 
and average over i, whereas we average over a and % to account for signals at both /* 
and 2/*; they follow Eq. (55) of Thorne (1987), who, in the context of bursting rather 
than continuous- wave sources, multiplies h c by (2/3) 1 ^ 2 to reflect a statistical preference for 
sources with directions and polarizations that give larger SNRs (because they can be seen 
out to greater distances); and they assume f c = 2/* instead of f c = 9/*/5 as required by 
(10). 
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4.2. Oscillations versus static mountain 

We now compare a star with an oscillating mountain against a star whose mountain is 
in equilibrium. We compute (10) and (11) directly from e(t) as generated by ZEUS-3D (see 
§2 and 3.1), i.e. without assuming that h + (f) and h x (f) are pure 5 functions at f — /*, 2f*. 

Table 1 lists the SNR and associated characteristic quantities for three M a values (and 
b = 10) for both the static and oscillating mountains. The case of a particular a and % 
(a — i — 7r/3) is shown along with the average over a and i (Thorne 1987; Bildsten 1998; 
Brady et al. 1998). We see that the oscillations increase the SNR by up to ~ 15 per cent; 
the peaks at / = /*, 2/* are the same amplitude as for a static mountain, but additional 
signal is contained in the sidebands. At least one peak exceeds the LIGO II noise curve in 
Figure 3 in each polarization. 



4.3. Detectability versus M a 

The SNR increases with M a , primarily because e increases. The effect of the oscillations 
is more complicated: although the Alfven sidebands increase in amplitude as M a increases, 
their frequency displacement from f — f* and / = 2/* decreases, as discussed in §3.2, so 
that the extra power is confined in a narrower range of /. However, e and hence the SNR 
plateau when M a increases above M c (see §3.1). The net result is that increasing M a by a 
factor of 10 raises the SNR by less than a factor of two. The SNR saturates at ~ 3.5 when 
averaged over a and % (multiple sources), but can reach ~ 6 for a particular source whose 
orientation is favorable. For our parameters, an accreting neutron star typically becomes 
detectable with LIGO II once it has accreted M a > 0.1M C . The base of the mountain may 
be at a depth where the ions are crystallized, but an analysis of the crystallization properties 
is beyond the scope of this paper. 
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Table 1: Signal-to-noise ratio 



/* [kHz] 


M a /1O" 4 M 


/ c [kHz] 


h c /io- 25 


SNR 




Static 


a = 7r/3 


% = 7r/3 




0.6 


0.16 


1.003 


0.83 


2.22 


0.6 


0.8 


1.003 


1.24 


3.34 


0.6 


1.6 


1.003 


1.35 


3.61 




Static 


0. 






0.6 


0.16 


1.08 


0.89 


2.22 


0.6 


0.8 


1.08 


1.33 


3.34 


0.6 


1.6 


1.08 


1.44 


3.61 




Oscillating 


a = 7r/3 


i = 7r/3 




0.6 


0.16 


1.008 


1.40 


2.63 


0.6 


0.8 


1.003 


2.15 


4.02 


0.6 


1.6 


1.004 


2.27 


4.25 




Oscillating 


0. 






0.6 


0.16 


1.056 


1.40 


2.45 


0.6 


0.8 


1.048 


2.14 


3.74 


0.6 


1.6 


1.048 


2.26 


3.95 
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5. DISCUSSION 

A magnetically confined mountain forms at the magnetic poles of an accreting neutron 
star during the process of magnetic burial. The mountain, which is generally offset from the 
spin axis, generates gravitational waves at /* and 2/*. Sidebands in the gravitational- wave 
spectrum appear around /* and 2/* due to global MHD oscillations of the mountain which 
may be excited by stochastic variations in accretion rate (e.g. disk instability) or magnetic 
footpoint motions (e.g. starquake). The spectral peaks at /* and 2/* are broadened, with 
full-widths half-maximum Af « 0.2 kHz. We find that the SNR increases as a result of 
these oscillations by up to 15 per cent due to additional signal from around the peaks. 

Our results suggest that sources such as SAX J1808.4— 3658 may be detectable by next 
generation long-baseline interferometers like LIGO II. Note that, for a neutron star accreting 
matter at the rate M a « 10 _11 M o yr- 1 (like SAX J1808.4-3658), it takes only 10 7 yr to 
reach S/N > 3. 3 The characteristic wave strain /i c ~4x 10~ 25 is also comparable to that 
invoked by Bildsten (1998) to explain the observed range of /* in low-mass X-ray binaries. 
An observationally testable scaling between h c and the magnetic dipole moment has 
been predicted (Melatos & Payne 2005). 

The analysis in §3 and §4 applies to a biaxial star whose principal axis of inertia coincides 
with the magnetic axis of symmetry and is therefore inclined with respect to the angular 
momentum axis J in general (for a^O). Such a star precesses (Cutler & Jones 2001), a fact 
neglected in our analysis up to this point in order to maintain consistency with Bonazzola 
& Gourgoulhon (1996). The latter authors explicitly disregarded precession, arguing that 
most of the stellar interior is a fluid (crystalline crust < 0.02M,), so that the precession 
frequency is reduced by ~ 10 5 relative to a rigid star (Pines & Shaham 1974). Equations 
(4) and (5) display this clearly. They are structurally identical to the equations in both 
Bonazzola & Gourgoulhon (1996) and Zimmermann & Szedenits (1979), but these papers 
solve different physical problems. In Zimmermann & Szedenits (1979), Q differs from the 
pulsar spin frequency by the body-frame precession frequency, as expected for a precessing, 
rigid, Newtonian star, whereas in Bonazzola & Gourgoulhon (1996), Q exactly equals the 
pulsar spin frequency, as expected for a (magnetically) distorted (but nonprecessing) fluid 
star. Moreover, 9 (which replaces a) in Zimmermann & Szedenits (1979) is the angle between 
the angular momentum vector J (fixed in inertial space) and the principal axis of inertia 



3 On the other hand, EOS 0748—676, whose accretion rate is estimated to be at least ten times greater, 
at M a > lO lo M0yr _1 , has /* = 45 Hz (from burst oscillations) and does not pulsate, perhaps because 
hydromagnetic spreading has already proceeded further (fj, < 5 x 10 27 Gcm~ 3 (Villarreal & Strohmayer 
2004). 
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e 3 , whereas a in Bonazzola & Gourgoulhon (1996) is the angle between the rotation axis 
O and axis of symmetry \x of the (magnetic) distortion. Both interpretations match on 
time-scales that are short compared to the free precession time-scale r p ~ (/*e) _1 , but the 
quadrupole moments computed in this paper (e ~ 10" 7 ) and invoked by Bildsten (1998) 
to explain the spin frequencies of low-mass X-ray binaries (10~ 8 < e < 10~ 7 ) predict r p of 
order hours to days. The effect is therefore likely to be observable, unless internal damping 
proceeds rapidly. Best estimates (Jones & Andersson 2002) of the dissipation time-scale 
give « 3.2yr(Q/10 4 )(0.1kHz//;) (J /10 44 g cm 2 ) (10 38 g cm 2 // d ), where J d is the piece of the 
moment of inertia that "follows" e 3 (not f2), and 400 < Q < 10 4 is the quality factor of 
the internal damping [e.g. from electrons scattering off superfluid vortices (Alpar & Sauls 
1988)]. 4 



4 Precession has been detected in the isolated radio pulsar PSR B1828— 11 (Stairs et al. 2000; Link & 
Epstein 2001). Ambiguous evidence also exists for long-period (<~ days) precession in the Crab (Lyne et al. 
1988), Vela (Deshpande & McCulloch 1996), and PSR B1642-03 (Shabanova et al. 2001). Of greater 
relevance here, it may be that Her X-l precesses (e.g. Shakura et al. 1998). This object is an accreting 
neutron star whose precession may be continuously driven. 
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Table 2: Precession scenarios and associated gravitational wave signals 





biaxial, e 3 | Q 


triaxial, e 3 | ft 




e 3|| A* 


zero GW 


GW at 2f* 


GW near /* and 2/* 




no precession 


no precession 


precession 




no pulses 


no pulses 


pulses 


e 3 It A* 


zero GW 


GW at 2/* 


GW near /* and 2/* 




no precession 


no precession 


precession 




pulses 


pulses 


pulses 



Note. — Here, e3 is the principal axis of inertia, fi is the axis of the magnetic dipolc, Jl is the spin axis, 
and /* is the spin frequency. Entries containing /* and/or 2/* indicate gravitational wave emission at (or 
near, in the case of precession) those frequencies; entries labelled 'zero GW indicate no gravitational wave 
emission. We also specify whether or not each scenario admits X-ray pulsations. 
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Some possible precession scenarios are summarized in Table 2. If we attribute persis- 
tent X-ray pulsations to magnetic funnelling onto a polar hot spot, or to a magnetically 
anisotropic atmospheric opacity, then the angle between /j, and fi must be large, leading to 
precession with a large wobble angle, which would presumably be damped on short time- 
scales unless it is driven (cf. Chandler wobble). Such a pulsar emits gravitational waves 
at a frequency near /* (offset by the body- frame precession frequency) and 2/*. However, 
the relative orientations of /x, fl, and e 3 are determined when the crust of the newly born 
neutron star crystallizes after birth and subsequently by accretion torques. This is discussed 
in detail by Melatos (2000). If viscous dissipation in the fluid star forces f2 to align with 
IX before crystallization, and if the symmetry axis of the crust when it crystallizes is along 
f2, then e3 (of the crystalline crust plus the subsequently accreted mountain), fi, and Q are 
all parallel and there is no precession (nor, indeed, pulsation). But if the crust crystallizes 
before fl has time to align with fi, then e 3 and Q are not necessarily aligned (depending 
on the relative size of the crystalline and pre-accretion magnetic deformation) and the star 
does precess. Moreover, this conclusion does not change when a mountain is subsequently 
accreted along fi; the new e 3 (nearly, but not exactly, parallel to /x) is still misaligned with 
fl in general. Gravitational waves are emitted at /* and 2/*. Of course, internal dissipa- 
tion after crystallization (and, indeed, during accretion) may force O to align with e 3 (cf. 
Earth). 5 ' 6 If this occurs, the precession stops and the gravitational wave signal at /* disap- 
pears. The smaller signal at 2/* persists if the star is triaxial (almost certainly true for any 
realistic magnetic mountain, even though we do not calculate the triaxiality explicitly in this 
paper) but disappears if the star is biaxial (which is unlikely). To compute the polarization 
waveforms with precession included, one may employ the small-wobble-angle expansion for 
a nearly spherical star derived by Zimmermann (1980) and extended to quadratic order by 
Van Den Broeck (2005). This calculation lies outside the scope of this paper but constitutes 
important future work. 

Recent coherent, multi-interferometer searches for continuous gravitational waves from 
nonaxisymmetric pulsars appear to have focused on the signal at 2/*, to the exclusion of the 
signal at /*. Examples include the SI science run of the LIGO and GEO 600 detectors, which 
was used to place an upper limit e < 2.9 x 10~ 4 on the ellipticity of the radio millisecond 
pulsar J1939+2134 (The LIGO Scientific Collaboration: B. Abbott et al. 2004b), and the 
S2 science run of the three LIGO I detectors (two 4-km arms and one 2-km arm), which 
was used to place upper limits on e for 28 isolated pulsars with /* > 25 Hz (The LIGO 



5 Accreting millisecond pulsars like SAX J1808.4— 3658 do not show evidence of precession in their pulse 
shapes, but it is not clear how stringent the limits are (Galloway, private communication). 

6 We do not consider the magnetospheric accretion torque here (Lai 1999). 
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Scientific Collaboration: B. Abbott et al. 2004a). Our results indicate that these (time- 
and frequency-domain) search strategies must be revised to include the signal at /* (if the 
mountain is static) and even to collect signal within a bandwidth Af centered at /* and 2/* 
(if the mountain oscillates). This remains true under several of the evolutionary scenarios 
outlined above when precession is included, depending on the (unknown) competitive balance 
between driving and damping. 

The analysis in this paper disregards the fact that LIGO II will be tunable. It is 
important to redo the SNR calculations with realistic tunable noise curves, to investigate 
whether the likelihood of detection is maximized by observing near /* or 2/*. We also do not 
consider several physical processes that affect magnetic burial, such as sinking of accreted 
material, Ohmic dissipation, or Hall currents; their importance is estimated roughly by 
Melatos & Payne (2005). Finally, Doppler shifts due to the Earth's orbit and rotation (e.g. 
Bonazzola & Gourgoulhon 1996) are neglected, as are slow secular drifts in sensitivity during 
a coherent integration. 

This research was supported by an Australian Postgraduate Award. 
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